A nonlinear Ramsey interferometer operating beyond the Heisenberg limit 
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We show that a dynamicaUy evolving two-mode Bose-Einstein condensate (TBEC) with an adi- 
abatic, time-varying Raman couphng maps exactly onto a nonlinear Ramsey interferometer that 
includes a nonlinear medium. Assuming a realistic quantum state for the TBEC, namely the SU(2) 
coherent spin state, we find that the measurement uncertainty of the "path-difference" phase shift 
scales as the standard quantum limit (l/\/iV) where A'^ is the number of atoms, while that for the 
interatomic scattering strength scales as , overcoming the Heisenberg limit of 1/N . 

PACS numbers: 03.75.Dg,03.75.Mn,03.75.Gg 



High-precision interferometry is one of the most im- 
portant tools of metrology enabling one to infer various 
properties of the system under consideration through the 
measurement of the phase shift. Ramsey interferometry 
provides a way to extract information about the changes 
in the system Hamiltonian H at time t via the phase 
shift (j) = H{t')dt' /h. In quantum systems, it is pos- 
sible to achieve measurement uncertainties approaching 
the Heisenberg limit A0 ~ where N is the number 
of particles conjugate to the phase variable (j) provided 
one uses carefully chosen entangled input states such as 
Schrodinger's Cat state. 

It has long been accepted that the Heisenberg limit is 
the ultimate limit to measurement; however recently it 
was shown using parameter estimation theory that mea- 
surement uncertainty of the order where k is the 
number of parameter-sensitive terms, is possible[l|. In 
particular, Bose-Einstein condensates (BECs) with two- 
body collisions may be able to achieve up to A(/) ~ 
accuracy in measurement of atom-atom interac- 
tions through a modulation the scattering length using 
Feshbach resonance or density variation due to gravita- 
tional gradients. In a related recent work[3|, the quantum 
limit to the measurement of atomic scattering length was 
studied by considering the Heisenberg exclusion principle 
applied to a squeezing Hamiltonian, and finding the opti- 
mal spin squeezed state generated using a separate time 
dependent Hamiltonian [3[. 

In this paper, we show that the temporal evolution of 
a TBEC such as Na atoms in the |F = l^Mp = ±1) 
hyperfine states trapped in an optical dipole trap with 
Raman coupling maps directly onto a nonlinear Ram- 
sey interferometer. A nonlinear interferometer, as op- 
posed to a normal (linear) interferometer, includes non- 
linear medium in one or both arms. Such TBEC systems 
have already been realized experimentally 0, Q, and as 
we show, can achieve measurement accuracy better than 
the Heisenberg limit. It is noted that nonlinear inter- 
ferometers have been studied previously (1, 0|, although 
never in the context of BEC. 

A quantum interferometer can be described in terms 
of the angular momentum operators as a transformation 
operator: 



The 50/50 beam splitter and the phase shifter are given 
by B± = exp{±iTTJx/2) and 'P((/>) = exTp{i(j)Jz) where 

Jx = ^{J+ + J-)i Jy = ^ J-)^ "^+(-) — ^1(2)^2(1) 

and Jz = ^{a\ai — a\a2) with ai and 02 being the two an- 
nihilation operators for the two input modes into interfer- 
ometer. For the TBEC with a Raman coupling like that 
considered here, the two annihilation operators ai and 
0,2 correspond to the atoms in the two hyperfine states. 
The time evolution operator U (t) for this system defined 

by m))=Uit)\m) isii: 

U{t) = R^e~'"''R, (2) 

where R = (^-^^--h)!-^ and H' = 2f7J^ - \Jl. Q is the 
tunneling coupling and q is the strength of the scattering 
interaction between the bosons. As shown earlier 0, the 
detuning of the laser from the transition between the 
two species is set to be zero to make the Hamiltonian 
diagonal in the Jz representation. This also prevents the 
generation of an additional geometric phase on top of the 
dynamical phase. 

The overall action of the time evolution operator U{t), 
Eq. ^ , can clearly be mapped onto a nonlinear Ramsey 
interferometer with the transformation operator 



(3) 



B- — W = exp(— i7rJy/2) and B+ = R are the 
two 50/50 beam splitters, while V{(j)'i) = e-*'^'i'^'- and 
^{4>'2) — e"^2'^j/2 represent, respectively, the "path- 
difference" phase shifter and the nonlinear medium. The 
phase shifts 0^ = 2Qt and (^2 = 1^ ^.re given by the Hamil- 
tonian dynamics. The presence of the intrinsic temporal 
evolution has to be taken into account in the measure- 
ment of (f)'i and 02 i-^- ^-^^y measurement will be shifted 
at the rate of 2f2 and q per unit time. Since we are inter- 
ested in measuring the changes in the phases (p^, i = 1,2, 
we shall redefine 0^ to be (p'l — 2Ut + (pi and 02 — qt + (t)2 
and concern ourselves with the measurement of pi. 

For our input state we shall consider an SU{2) atomic 
coherent state or a coherent spin state (CSS), \6, p)) which 
is a reasonable quantum state representing a TBECfOj, 



(1) p. It is noted that R{t)\6, 0) = ^m(^ + \ 



2 



At)\j,m) where 72.^^(6', 0) is defined 
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(4) 



Since the azimuthal angle (p simply to shifts the origin, we 
shall only consider CSS with (j) = Q in this paper. Exotic 
input states such as the NOON or the Yurke stateQ 
will be considered elsewhere as they are currently not 
yet practical in the context of TBEC. 

The simplest possible scenario is to measure the "path- 
difference" phase shift (j)i while applying a magnetic field 



to tune 



via the Feshbach resonance i.e. no nonlin- 



ear perturbations to the Hamiltonian. This is the stan- 
dard Ramsey interferometry which has been studied ex- 
tensively. The fact that a TBEC is used instead of the 
thermal atoms simply provides clean signals owing to the 
inherent long range coherence of a condensate. We will 
not consider this case any further in this paper. What's 
more interesting is the case of finite q. Here the presence 
of the nonlinear component modifies the interferometric 
outcome 0i , and brings to the fore the question of the un- 
certainty associated with measuring the scattering length 
or (/)2. 

First, we analyze TBEC as a nonlinear Ramsey inter- 
ferometer in the idealized situation where the measure- 
ment of the phase is carried out as a projective mea- 
surement onto a phase state. We estimate the measure- 
ment uncertainty using the Cramers-Rao inequality in 
such cases. Then a more practical scheme, measurement 
of the atom number difference as a function of the phase 
shifts 01 and (j)2 , is considered along with the correspond- 
ing measurement uncertainties. The fundamental limit 
to the phase shift measurements can be calculated by 
first defining the positive valued operator measure E{(j)) 
such that the probability density of the corresponding 
measurement result is: P{4>) = Tr[/5i?((/))] where p is the 
density matrix for the system. As in Ref. pT| . we define 
the normalized phase state 

|j,<i>) = (2j- + l)-i/2 ^ e™^*|j,m.). (5) 

SO that E{<^)d^ = (2j + 1)| j, $) (j, $|d$/27r, and for an 
arbitrary input state the probability density of the 
measurement result is 



(6) 



With a CSS input, the phase measurement gives a 
probability density distribution 
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FIG. 1; Probability density P($) for the initial CSS \e = 
7r/4, = 0) (Solid line) and \e = 7r/2, = 0) (Dashed line) at 
different times fit = O.TSvr/, 2.57r, 60.257r (Top, middle, bot- 
tom rows respectively). Left column: q = 0; right column 
g/0. 



where 
matrix: 



(7r/2) =^ {j,rnz\j,mx).j; is the Wigner d- 



(j - m^)!(j + TO2)! 



-,1/2 



(7) 



for TUz — urix > — 1 and mz + rux > —1. P'--°''^\x) de- 
note the Jacobi polynomials. Symmetries give „ = 

(-ir^7'"^rf^™.,™„ We plot the probability 

density in Fig. [T]at various times starting from the initial 
states 16* = 7r/4) (Sohd line) and \e (Dashed line). 

The initial Dicke state \6 = 0) is not considered as it is 
orthogonal to the projective measurement on the phase. 
In order to highlight the effect of nonlinearity on the 
measurement of 0i , we plot in the left column the case of 
q = for comparison with the corresponding results with 
q 7^ in the right column. We choose q = 3/N which 
corresponds to the Josephson regime^]. The presence of 
nonlinearity generally degrades the performance of the 
interferometer as evidenced by the increase in the width 
of the probability distribution. It is also notable that the 
probability density becomes multiple peaked after a long 
time. This may be interpreted as the generation of a su- 
perposition state due to nonlinearity as studied by Yurke 
and StolerjT^. It is clear that to use the TBEC as an 
effective interferometer based on projective measurement 
onto phase states, q needs to be minimzed and the time 
of measurement must be kept relatively short. 

The uncertainty in phase measurement can be stud- 
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FIG. 2: Time evolution of logjv[A<l>] for the initial CSS 
\e = TT/4,(f> = 0) (Solid line) and \e = TT/2,(f) = 0) (Dashed 
line) . Top row: direct calculation from the probability density 
P("I>); Bottom row: Cramers- Rao lower bound. Left column: 
q = 0; right column q ^ 0. 
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FIG. 3: { Jz) jN as a function of the changes in the phase 
shifts (^1 (Left column) and 02 (Right column) for the initial 
CSS |6I = 0, = 0) (Dashed line) and = 7r/4,<^ = 0) (Solid 
line) at different times f2t = 7r/4, tt, Gtt (Top, middle, bottom 
rows respectively). 



ied using the standard techniques of probability theory, 
particularly the Cramers-Rao lower bound (CRLB). The 
CRLB establishes the lower bound on the phase shift esti- 
mate where the phase uncertainty scales as A<J> — \j \/T\x 
where is the Fisher information defined by 



1 

2^ 



lnP($) 



P($)d$. 



(8) 



In Fig. [21 we plot the quantity logjY[A$] where iV is 
the total number of atoms and A$ is the uncertainty in 
phase. We used the standard deviation A<I> calculated 
directly from the probability distribution P(<i>) (Fig. [T]) 
and the CRLB, where the CRLB effectively gives a time 
averaged value of the directly calculated uncertainty. It 
is noted that in all these figures A$ > l/A^^^^, the stan- 
dard quantum limit. It is also noted that although the 
uncertainty associated with the 16* — 7r/2) state is lower 
than that of the \d = 7r/4) state ior q — it quickly 
loses this advantage with q > 0, indicating sensitivity to 
dephasing due to interatomic collisions. The 16* ~ 7r/4) 
state is therefore a more robust state for interferometry 
in the presence of nonlinearity. 

Next, instead of projective measurement onto a 
phase state, we consider projective measurement of 
the atom number difference. The total number of 
atoms measured indicates the number of "input" atoms 
while the atom number difference, {Jz{4''ii4''2)) = 
(i/'(O)|i:t(0;,0^)j^i:(^;,^^)|^(O)), allows us to infer the 
phase shift and is equivalent to measuring the number 
of atoms at each of the output port of a typical Mach- 
Zehnder interferometer. An analytic expression for {Jz) 



is given byQ: 

(>/.K,02)> 

where we 



N/2~l 

^ X'(6l,m)tan"^ 

m=-N/2 



e-TT/2 



X COS 



(9) 



/-<Af ( N_ 

^JV/2+m+l V 2 



have defined 

71+1) COS^^ 



Figure [3] shows ( Jz) /N as a function of the changes in the 
phase shifts (pi and (j)2 for the initial states 16* = 0, </) = 0) 
and 16* = 7i'/4) at different times Vlt = 'k/A,'k,Q'k. In 
contrast to the earlier phase state projection method, 
the initial CSS 16* = 7i'/2) is known as a "self-trapping" 
state in the new context of projective number measure- 
ment, and gives trivial results. Since the interferometry 
is carried out at fixed times, we see in Fig. [3] clear 
sinusoidal fringes, without the "collapses and revivals" 
typical of temporal evolution. Even when measuring 
(j)2 we see clear fringes for a range of values around 
— 3g . . . Bq for the initial state \9 = 7r/4). This is possible 
because, for this choice of 6*, the factor ViO = 7r/4,TO) is 
narrow enough to limit the interfering effect of summing 
up the cosine terms. On the other hand, ViO = 0, m) is 
wider and the resulting interference fringes do not allow 
for a sensitive detection of small variations in 02 • 

Finally, we consider the phase resolution for this 
scheme which is given by: 



[AJ, 



|9(j.)/a0fc| 



A: = 1,2 



(10) 
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FIG. 4; Top: log^yiA^fc] plotted as a function of the angle 9 
of the initial CSS with iV = 1000. Solid line: k = 1, Dashed 
line: k — 2. Bottom: same quantity plotted as a function 
of the number of atoms, A'^. Solid line: k = 1 with 9 = 0. 
Dashed line: k = 2 with 9 = tt/4. 
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[oj , the variance 

with (J2(t)) = 

7r/4)(iV/2 - m) + 



where, as found in 
is [Air = {Jl) 
lEltN/2mm)t^n 

V{e,m){N/2 + m + l)+Y."n= 
7r/4)(iV/2 - TO - 1)] cos[2(j)\ - 2(t)'^{m + 1)]. Since the 
denominator of Eq. (|10p involves a function of the form 
sin — (j)'2 (to + ^)] , the quantity [A(/)fc]^ is minimized 
for the values of —(j)'2 {m + ^) — ±7r/2. This indicates 
that the measurement accuracy is dependent on the 
measurement values, where results such as = ±n/2 
and 02 = give optimum results. In particular, for 
a large number of atoms N one can approximate the 
coefficient V{e,m) by v^TV/Tre-f^^-JVsine)^^ and 
replace the sums by integrals / x)dx — N and 



jNV{9,x)dx = JxV{e,x)dx = N"^ . This leads to 
[AJ^]2 - aiV + /37V2 and |9(J^)/9<^i|2 ~ 77V2 
\d{Jz) I d(j)2\^ ^ ^N'^ where a,P,j are constants so that: 



[aN 



-2k+l 



-2/0+2] 



(11) 



where fc = 1,2. For k — 1 one has A0i ~ i.e. 
the standard quantum limit in accuracy for the measure- 
ment of (f>i. On the other hand, it is remarkable that 
with k = 2 i.e. measurement of the phase shift due to 
the interatomic interactions, A02 ~ l/N^^"^ < 1/N, im- 
plyingthat, although not reaching the theoretical limit of 
1/A^2[l|, such measurement for this CSS input state has 
uncertainty below the Heisenberg limit. We have verified 
this estimate numerically; the quantity logjy Acjjk calcu- 
lated as a function of the initial angle of the CSS, 9 at 
the optimal values of and 02 is plotted in Fig. SI The 
solid and the dashed line represent the uncertainty in the 
measurement of 0i and 02 respectively. It is clear that 
the best result is obtained for 6* = for the measurement 
of 01 and 9 = tt/A for the measurement of 02. In the 
bottom panel, we plot the result as a function of atom 
numbers for these chosen values of 9. It shows that the 
result is independent of number of atoms and, on aver- 
age, A01 l/TV^/^ and A02 \/N'^/^ which is indeed 
very close to the above estimate. 

In summary, we have shown that a TBEC with a 
Josephson-like coupling directly maps onto a nonlinear 
Ramsey interferometer where the phase shifts due to 
linear and nonlinear variations in the Hamiltonian are 
measured. The system is already experimentally avail- 
able and the state we consider is the realistic coherent 
spin state rather than some exotic quantum state. It 
was found that projective phase measurement reaches 
the standard quantum limit in accuracy while, remark- 
ably, projective number measurement of the phase shifts 
due to interatomic interactions were found able to over- 
come the Heisenberg limit, suggesting new implications 
for quantum metrology. 
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